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Abstract 

Empirical  evidence  shows  that  the  rate  of  zero-crossings  of  many 
stochastic  processes  tends  to  increase  by  repeated  differencing.  This 
motivates  the  definition  of  a  class  of  processes  whose  expected  oscilla¬ 
tion  increases  monotonical ly  by  repeated  differencing.  The  class  of 
strictly  stationary  processes  is  a  subclass  of  this  class.  It  is  shown 
that  there  is  a  limit  to  oscillation  by  proving  that  the  point  processes 
of  zero-crossings  obtained  by  repeated  differencing  converge. 
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STRICTLY  OSCILLATORY  PROCESSES 


by 

B.  Kedem  and  D.  Martin 


1 .  Introduction 

In  this  paper  we  introduce  a  class  of  random  processes  to  which  we 
refer  as  strictly  oscillatory  and  suggest  a  method  to  monitor  the  oscilla¬ 
tion  observed  in  such  processes.  When  a  process  is  second-order  sta¬ 
tionary,  the  oscillation  observed  in  the  process  is  described  very 
effectively  by  the  spectrum.  When  the  process  is  nonstationary,  various 
attempts  have  been  made  to  extend  the  notion  of  the  spectrum  to  model 
the  time  varying  spectral  content  of  the  process.  However,  a  process 
need  not  possess  moments  at  all  and  still  appear  to  be  oscillatory. 

What  is  needed  then  is  a  way  to  describe  oscillation  in  random  phenome¬ 
na  removed  from  stationarity  assumptions  and  independent  of  any  moment 
conditions. 

In  many  respects  the  simplest  way  to  describe  the  oscillation  ob¬ 
served  in  a  stochastic  process,  stationary  or  nonstationary,  is  through 
the  point  processes  obtained  from  higher  order  crossings.  The  advan¬ 
tages  offered  by  such  zero-crossing  counts  are  as  follows: 

1.  The  pattern  of  oscillation  changes  can  be  detected  and  described 
directly  by  zero-crossing  counts  without  recourse  to  any  Fourier 
analysis.  Thus  we  gain  simplicity. 

2.  The  zero-crossing  counts  observed  in  finite  series  in  discrete  time 
possess  all  moments  regardless  of  whether  the  original  process  has 
moments  or  not.  Thus  we  relax  the  requirement  of  finite  moments. 
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3.  In  scientific  and  engineering  applications  the  amounts  of  data  can 
be  very  large  and  at  times  forbidding.  Zero-crossing  counts  of 
filtered  series  provide  effective  data  reduction. 

4.  As  we  shall  see,  it  is  possible  to  construct  useful  graphical  plots 
from  the  rates  of  higher  order  crossings  to  monitor  spectral  changes. 
In  this  regard  the  monotonicity  of  higher  order  crossings  plays  a 
crucial  role. 

We  shall  show  that  the  class  of  strictly  oscillatory  processes  is 
quite  large,  and  prove  a  certain  weak  convergence  associated  with  such 
processes.  This  convergence  means  that  in  some  sense  there  is  a  limit 
to  oscillation. 

Sections  2  and  3  motivate  our  definition  of  strictly  oscillatory 
processes  given  in  section  4.  In  section  5  we  prove  the  convergence  of 
point  processes  defined  by  higher  order  crossings,  while  section  G  dis¬ 
cusses  briefly  the  Gaussian  case. 
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2.  Evolutionary  Spectra  Approach  to  Oscillatory  Processes 

An  attempt  to  model  oscillatory  processes  has  been  made  in  Priestley 
(1981,  ch.  11).  It  is  important  to  make  a  distinction  between  the  two 
approaches  and  to  point  to  the  difficulties  one  encounters  in  trying  to 
extend  the  theory  of  second-order  stationary  processes  to  the  nonsta¬ 
tionary  case.  It  should  be  borne  in  mind  that  our  approach  is  completely 
different  and  is  graphical  combinatorial  in  nature. 

In  classical  spectral  analysis  of  time  series,  it  is  generally 
assumed  that  time  series  are  generated  by  second-order  stationary  pro¬ 
cesses  where  the  autocovariance  function  and  the  spectrum  form  a  Fourier 
pair.  When  stationarity  breaks  down,  this  Fourier  relation  is  no  longer 
true  since  the  autocovariance  is  time  dependent.  There  were  several 
attempts  to  define  time-varying  spectra  of  nonstationary  processes. 
Perhaps  the  most  complete  approach  is  that  of  "evolutionary  spectra" 
developed  by  Priestley  (1981).  The  basic  motivation  underlying  the 
evolutionary  spectra  approach  stems  from  the  need  to  model  the  local 
behavior  of  nonstationary  processes.  Whereas  the  spectrum  of  a  sta¬ 
tionary  process  describes  the  power- frequency  distribution  for  the 
whole  process,  the  evolutionary  spectrum  is  time  dependent  and  de¬ 
scribes  the  local  power- frequency  distribution  at  each  time  instant. 

A  stationary  process  { Z ^ }  with  spectral  distribution  function  F(ui) 
admits  the  spectral  representation 


eltu)dr(o)),  t  =  0,±1,..., 


(1) 


'-u 


where  £(w)  is  a  process  of  orthogonal  increments  and  satisfies 


E|d^U)|2  =  dF(w) . 


w.  /?,<•.  .r.y.V.  ■ 


V- 
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This  representation  allows  us  to  identify  the  component  in  which  has 
frequency  iu  and  to  determine  the  contribution  of  this  component  to  the 
total  power  of  the  process.  A  way  to  define  a  class  of  nonstationary 
processes  is  to  replace  the  exponential  function  in  (1)  by  a  more 
general  class  of  functions.  Priestley  (1981,  ch.  11)  defines  an  oscil¬ 
latory  process  to  be  a  process  which  admits  the  representation 

[TT  •  . 

Z.  =  e  WA  (u))d£(uj),  t  =  0,±1,... 

—  7T 

where  for  each  w,  At(oj)  has  a  generalized  Fourier  transform  whose  modu¬ 
lus  has  an  absolute  maximum  at  the  origin,  and  F,(gj)  is  as  in  (1).  The 
evolutionary  spectrum  at  time  t  is  defined  as 

2 

dH^(oj)  =  [ A^. { oj )  |  dF(w),  -tt  <  w  <  ti 

and  it  represents  the  spectral  content  of  the  process  at  time  t.  From 
(2)  and  (3)  it  is  possible  to  construct  a  theory  of  prediction  and  fil¬ 
tering  for  nonstationary  processes  which  parallels  that  for  stationary 
processes. 

Though  useful  in  some  respect,  this  theory  has  its  drawbacks  when 
viewed  from  the  point  of  view  of  applications.  First,  the  representa¬ 
tion  (2)  is  not  unique  and  one  is  faced  with  the  problem  of  choosing 
At(uj).  A  possible  remedy  is  suggested  in  Priestley  (1981)  when  At(ui) 
is  a  slowly  varying  function  of  t  for  each  fixed  w,  but  is  too  techni¬ 
cal  to  be  described  here.  Second,  it  is  difficult  to  prove  that  a 
process  encountered  in  practice  can  in  fact  be  represented  by  (2). 
Third,  there  are  processes  which  appear  to  be  very  oscillatory  but 
possess  no  moments  or  just  first  order  moments  but  for  which  (2)  is 
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clearly  unsuitable.  Fourth,  when  using  the  evolutionary  spectrum  (3) 
to  detect  spectral  changes  in  nonstationary  time  series,  we  need  to 
consider  and  compare  many  different  pictures  or  graphs,  one  for  each 
time  point  t.  Thus  graphical  detection  of  spectral  changes  can  be 
quite  cumbersome  using  evolutionary  spectra.  Also,  the  estimation  of 
numerous  spectra  can  be  very  time  consuming.  A  less  stringent  approach, 
but  one  which  still  characterizes  the  local  character  of  nonstationary 
processes,  can  alleviate  these  difficulties.  Such  an  approach  can  be 
developed  by  higher  order  crossings  of  nonstationary  processes  and  is 
evidently  not  tied  to  any  moments  assumption. 


r, 

r 


a 


I 


3.  Zero-crossing  Rate  Processes 

Let  (Zt),  t  =  0,±1,...,  be  a  stochastic  process  defined  on  some 
probability  space  (ft,A,p).  it  is  assumed  the  process  is  not  a  constant 
with  probability  one  but  that  it  may  be  stationary  or  nonstationary  and 
it  may  or  may  not  possess  moments  of  any  order.  Let  V  be  the  backward 
difference  operator,  vZt  =  Zt  -  Zt  l,  VJZfc  =  vfV^Zj.),  and  define 
a  sequence  of  binary  processes  {X^}  by 

1. 


Xtj)  5  < 


Vj_1Zt  >  0 


0,  otherwise,  t  =  0,±1 ,, 


The  symbol  changes  in  {xj.^}  correspond  to  zero-crossings  in  discrete 
time  in  {V^Z^.}.  We  are  interested  in  the  number  of  such  zero- 
crossings  in  finite  time  series.  Define 

d(j)  =  I[X(j)  ?  X[jj],  t  =  0 ,±1 , . . . , 


where  I  [ ■  ]  is  the  indicator  function.  When  d”^  =  1  we  say  that  a 
zero-crossing  occurs  at  time  t  in  f V J ~ ^ Z ^ } .  The  number  of  zero- 


crossings  in 


VJ_1Zt,  t  =  1,2,...,N 


is  denoted  by  Dj  N  and  is  therefore  given  by 


V  d 
t=l 


(j) 


j  =  1,2,...  . 


The  Dj  ^  are  called  the  (number  of)  higher  order  crossings  in 
t  =  1,2,...,N,  or  HOC  for  short.  Higher  order  crossings  were  intro¬ 


duced  by  Kedem  and  Slud  (1981,  1982)  for  the  purpose  of  discrimination. 


while  their  applications  in  spectrum  analysis  are  discussed  in  Kedem 
(1986a,  1986b,  1987).  A  review  of  the  theory  and  application  of  HOC 
can  be  found  in  Kedem  (1986c).  This  work  demonstrates  that  when  it 
comes  to  stationary  Gaussian  processes,  HOC  possess  a  surprising  amount 
of  spectral  information.  In  fact,  in  this  case  the  sequence  of  expected 
HOC  determines  uniquely  the  correlation  function  fp^}  and  the  normalized 
spectral  distribution  function  F(oj)  and  we  write 

{EDjiN}  <=>  (Pk}  <=>  F(u)).  (4) 

Relation  (4)  emboldens  us  to  believe  that  higher  order  crossings  may  be 
helpful  in  tracing  spectral  changes  in  nonstationary  processes  as  well. 

Observe  that  in  the  stationary  case  the  zero-crossing  rate  D-  W/N 
converges  to  a  constant  in  some  sense  for  fixed  j  as  N  +  On  the 
other  hand  for  a  nonstationary  process  no  such  convergence  is  expected 
since  the  spectral  content  or  oscillatory  pattern  of  such  a  process 
may  be  time-varying.  This  motivates  our  definition  of  rate  processes 
for  monitoring  spectral  changes  and  changes  in  oscillation  patterns. 

Definition  1 .  A  (HOC)  rate  process  of  order  j  is  defined  for  each 
fixed  j  by  the  process 

D.  N 

Dj(N)  =  N  =  1,2,...  . 

It  is  seen  that  for  each  fixed j »  D j ( N )  traces  the  zero-crossing  rate 
of  { V J “ 1 Z t } ,  t  =  1,...,N,  as  a  function  of  N.  Evidently 

0  <  D  j  ( N )  <  1 

and  so,  any  rate  process  possesses  all  moments.  It  is  helpful  to  look 
at  the  graphs  of  rate  processes  obtained  from  different  time  series. 


The  figures  below  feature  rate  processes  corresponding  to  sta¬ 
tionary  and  nonstationary  processes.  Each  figure  displays  the  first 
eight  rate  processes  D1 ,(•),••• .Dg( •) ,  where  D^.)  is  the  lowest 
curve  while  Dg(*)  is  the  highest.  It  is  seen  that  changes  in  the 
spectral  content  and  in  the  oscillatory  pattern  in  each  case  are  well 
captured  by  the  rate  processes.  Also,  the  rate  processes  display  a 
monotone  convergence  towards  a  limiting  rate  process  which  can  be  in¬ 
terpreted  to  mean  that  there  is  a  limit  to  the  oscillation  displayed 
by  a  process  stationary  or  not.  It  is  this  monotone  behavior  observed 
in  numerous  phenomena  which  prompted  this  investigation  and  conse¬ 
quently  our  definition  of  strictly  oscillatory  processes. 

We  close  this  section  by  noting  a  useful  feature  of  the  figures. 

As  mentioned  earlier,  if  one  tries  to  plot  evolutionary  spectra  on  a 
single  coordinate  system  the  graphs  may  overlap  heavily  so  that  separate 
plots,  up  to  one  graph  for  each  time  point,  are  necessary  in  general. 

On  the  other  hand  our  figures  provide  a  compact  way  to  trace  the 
oscillatory  history  of  a  process  on  a  single  frame  and  overlapping  is 
prevented  due  to  the  displayed  monotonicity  of  rate  processes.  Now, 
if  one  tries  to  imitate  the  idea  of  plotting  other  quantities  as  func¬ 
tions  of  time,  overlapping  graphs  may  occur.  For  example.  Figure  12 
shows  the  autocorrelations  of  the  differences  VJZt  as  functions  of 
time  where  the  time  series  comes  from  the  utterance  of  "six".  It  is 
seen  that  both,  the  rate  processes  and  the  autocorrelation  processes, 
trace  the  spectral  changes  in  the  time  series  but  in  the  rate  processes 
no  overlapping  occurs. 
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Figure  1.  A  time  series  from  the  utterance  of  "zero"  and  the  corre¬ 
sponding  first  eight  rate  processes. 


r«.4». 


l.l  1.4^ 


■nmni 


10 


‘■Pl.FXH  PATA  'JNF. 

I00P  da  i  a  poin  r  nur  of  0000 


-2000  — |— i ; - fl  ,  -| - p 

0  1 000 


ft  -j-  t  1  -r  r  1  T 
.’000  -001* 


T— f  I- [  7  r  -I  -f- 1  I  - T—  r  1 
4000  '1000  r,000 


l.)  I  1 


Figure  2.  A  time  series  from  the  utterance  of  "one"  and  the  corre 
sponding  first  eight  rate  processes. 
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A  time  series  from  the  utterance  of  "two"  and  the  corre 
sponding  first  eight  rate  processes. 
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A  time  series  from  the  utterance  of  "four"  and  the  corre 
sponding  first  eight  rate  processes. 
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A  time  series  from  the  utterance  of  "five"  and  the  corre 
sponding  first  eight  rate  processes. 
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Figure  11.  The  rate  processes  of  the  nonstationary  oscillatory  process 
(Me lard  (1985)), 
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Rate  processes  and  correlation  process  from  an  utterance 
of  "six".  In  the  correlation  processes  overlapping  occurs 


4.  Definition  of  Strictly  Oscillatory  Processes 

The  monotone  behavior  of  rate  processes  observed  in  the  figures  is 
quite  a  general  phenomenon  and  can  be  attributed  to  the  difference  opera¬ 
tion.  For  example,  in  the  second-order  stationary  case,  a  difference 
operation  acts  as  a  high  pass  filter  which  gives  more  power  to  higher 
frequencies.  This  amplification  causes  the  process  to  be  more  and  more 
oscillatory  with  each  application  of  V.  In  general,  differencing  causes 
values  at  neighboring  points  to  be  negatively  correlated  and  thus  in¬ 
crease  the  number  of  symbol  changes  in  {X^}.  It  seems  therefore  that 
the  requirement  that  ED-(N)  increases  as  j  -+  °°  for  each  fixed  N  is 

J 

quite  reasonable  due  to  this  observed  generality.  This  leads  to  our 
definition  of  strictly  oscillatory  processes. 

To  define  the  class  of  strictly  oscillatory  processes,  and  at  the 
same  time  describe  the  monotone  convergence  observed  in  the  figures  we 
appeal  to  the  theory  of  point  processes  as  developed  by  Kallenberg  (1976). 

Recall  that  { Zt } ,  t  =  0,±1,...,  is  a  process  defined  on  the  proba¬ 
bility  space  (fl,A,P)  and  let  G  =  (0,°°)  be  equipped  with  the  Borel 
field  G.  Further  let  6  be  the  class  consisting  of  all  bounded  sets  in  G. 
The  set  of  all  locally  finite  measures  p  that  are  nonnegative  and  integer 
valued  is  denoted  by  M.  Thus  p  c  N  if  for  B  e  t> 

p(B)  e  Z+  ~  {0,1,2,...}. 

In  W  we  introduce  the  a-algebra  F  generated  by  the  mappings  p  -+  p(B), 

B  e  6.  This  is  the  smallest  o-algebra  that  makes  these  mappings  measur¬ 
able.  By  a  point  process  on  S  we  mean  any  measurable  mapping  of  ($i,A,P) 
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With  an  obvious  extension  of  the  previous  notation,  let  D.(B), 

J 

B  e  8,  be  the  number  of  zero-crossings  generated  by  in  B.  Then 

for  each  realization  of  { } ,  Dj(*)  defines  a  measure  in  N  and  so  for 
each  fixed  j,  Dj*)  is  a  point  process.  In  other  words,  the  values  of 

v) 

D-(*)  when  ranging  over  fi  are  measures  in  W.  The  sequence  can 

J  J 

now  be  used  in  defining  the  class  of  oscillatory  processes. 

The  following  definition  reflects  our  experience  that  by  differencing 
any  finite  section  of  a  process  the  corresponding  higher  order  crossings 
tend  to  increase.  Note  that  since  t  is  defined  over  all  the  integers, 
theoretically  we  do  not  have  any  problem  with  repeated  differencing  ap¬ 
plied  to  finite  sections.  On  the  other  hand,  in  practice  only  finite 
records  are  available  and  we  lose  one  observation  with  each  difference. 
However  this  does  not  pose  any  serious  problem  since  the  amount  of  dif¬ 
ferencing  needed  in  practice  is  very  small  as  is  well  indicated  by  the 
above  figures. 

Definition  2.  A  stochastic  process  { Z t } ,  t  =  0,tl,...,  is  called 
strictly  oscillatory  if  and  only  if 

EDj (B)  <  EDj+1(B) 

for  every  B  e  8  and  j  =  1,2 . 

The  class  of  strictly  oscillatory  processes  is  quite  large. 

Strictly  stationary  processes  are  members  of  this  class  and  second-order 
stationary  processes  become  strictly  oscillatory  by  undergoing  a  sufficient 
amount  of  differencing.  However,  oscillatory  processes  in  the  sense  of 
Priestley  (1981)  need  not  be  strictly  oscillatory.  Our  first  example  con¬ 
cerns  stationary  Gaussian  processes. 
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Theorem  1 .  A  zero  mean  stationary  Gaussian  process  is  strictly  oscilla¬ 
tory. 

Proof.  Let  p.j ,  Py(l)  be  the  first  autocorrelations  in  Z  and  VZ  respec¬ 
tively.  Then  from  Kedem  (1984)  PyO)  <  p^ .  By  stationarity,  it  is 
sufficient  to  consider  the  sequence  of  higher  order  crossings.  Consider 
D-|  N  and  D£  Then  the  Gaussian  assumption  entails 

ED1,N  =  N4  -  7  sin"1p1)  <  N(^  -  1  sin_1Py(l))  =  ED^. 

But  since  the  differences  { Zt }  are  again  stationary  and  Gaussian  with 
mean  zero  we  also  have  ED,,  N  <  ED^  and  ED^  N  <  ED^  ^ ,  etc.  There¬ 
fore 

EDj,N  "  EDj+l ,N*  a11  J  =  1 ,2 . 

Since  the  distributions  in  { Z^ }  are  unaffected  by  time  shifts  due  to 
stationarity,  the  proof  is  complete.  a 

We  can  actually  prove  a  stronger  result  of  which  Theorem  1  is  only 
a  special  case.  This  pertains  to  strictly  stationary  processes  in 
general.  In  this  case  no  moments  assumption  is  needed  at  all. 

Theorem  2.  Every  strictly  stationary  process  {Z^.}  is  strictly  oscillatory. 
Proof.  Again,  by  stationarity  it  is  sufficient  to  consider  the  HOC  D-  M. 

J  ,  11 

Define 

xjj)  i  P(x[j)  =  l  j  =  1  ),  p(j)  -  P(x[j)  =  1  ). 

EDj>N  =  2Np(j)(l  -  xjj)). 


Then 


Now,  in  general  we  have  (See  Lemma  1  below) 


Dj  »N  5  °j+l,N  +  1 


surely.  Therefore  EDj  ^  <  ED j  +1  N  +  1  and  so 


2Np^(l-A^)  2Np^+1^(l  -  xjj+1^) 

—  <  - L_ —  +  1 


N 


N 


N 


As  N  «>  we  obtain  the  inequality 


Aj) 


2pVJ/(1  -  A  J/)  <  2p 


(j) 


,(j+l) 


(1-A 


(j+D- 


Therefore 

EDjfN  =  2Np(j)(l  -  a}j))  <  2Np(j  +  1)(l  -  AjJ  +  1))  =  EDj  +  ] ^  . 

Since  by  stationarity  every  stretch  of  data  of  fixed  length  is  endowed 
with  the  same  distribution,  it  follows  that  the  process  is  strictly 
osci 1 latory. 


Corollary  1.  A  sequence  of  independent  and  identically  distributed  ran 
dom  variables  is  strictly  oscillatory. 


Observe  that  the  original  sequence  may  be  completely  random  but 
that  the  higher  order  differences  become  highly  dependent  and  by  the 
corollary  also  more  and  more  oscillatory. 

Our  next  example  concerns  weakly  stationary  processes.  Here  we 
appeal  to  the  spectral  properties  of  such  processes  as  expressed  by 
their  spectral  measures.  But  first  we  prove  a  basic  lemma  which  builds 
on  the  interesting  relation  between  {xj^}  and  {X^  +  ^}. 
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Lemma  1 .  For  any  process  {Z^.}  there  exists  a  tQ,  1  <  tQ  <  N  and  tQ 
depends  on  j,  j  =  1,2,...,  such  that  for  every  N 


i  <  d< 

t=l  1  t=l  1 


(j+n 


with  probability  one. 

Proof.  Fix  j.  Suppose  Dj  ^  =  0  ;  then  we  are  done.  So  suppose 

Dj  N  >  0.  Then  there  is  at  least  one  symbol  change  in  {x|Jb,  t=l,...,N. 

Note  that  in  general  we  have  the  implication 

{X«j  t  x[j)}  c  {x[J  +  1)  =  x[j)}’  (6) 

It  follows  that  there  exists  t-,  such  that  xfJ+1^  =  xiJ'^.  But  (6)  implies 
that  each  symbol  change  in  {xj.J'}  produces  at  least  one  symbol  change  in 
(X^+1h  for  t-|  <  t  <  N.  On  the  other  hand,  again  by  (6)  and  due  to  an 
end  effect,  the  number  of  symbol  changes  in  (xj.J+1h  always  exceeds  that 
in  {x[j)}  minus  one  for  1  <  t  <  t-j .  It  follows  that  there  exists  a  tg 
which  depends  on  j  such  that  (5)  holds.  □ 

The  lemma  says  that  in  {xj.^+^},  1  <  t  <  N,  there  are  N-l  loca¬ 

tions  which  give  at  least  as  many  symbol  changes  as  the  corresponding 
locations  in  { X ^ J ^ } ,  while  no  information  is  available  about  the  re¬ 
maining  location  which  we  denote  by  tg.  Clearly  tg  may  change  with  j. 

Note  that  at  tg  we  may  lose  a  symbol  change  with  each  difference  as  is 
the  case  with  polynomials  in  t.  For  by  differencing  a  polynomial  of 
degree  n  we  obtain  another  polynomial  of  degree  n-l. 


V. 
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Theorem  3.  Let  { X^. }  be  a  weakly  stationary  process  with  spectral  dis¬ 
tribution  function  F.  Assume  that  it  is  in  the  spectral  support.  Then 
{ VJZt }  becomes  strictly  oscillatory  as  j 

Proof.  Again  by  stationarity,  it  is  sufficient  to  consider  the  higher 
order  crossings.  The  normalized  spectral  measure  of  { V J Z t }  is  given  by 

(sin  iu>)^dF(w) 

v  <(dw)  =  — - ; - 

f  (sinix)2JdF(X) 

-IT 

and  from  Kedem  and  Slud  (1982)  we  know  that  v,  converges  weakly  to  a 

J 

measure  supported  at  ±n.  That  is, 

V  •  =>  Jr  6  +  7T  6 

J  2  -TT  2  TT 

where  6  denotes  the  point  mass  at  u.  Therefore  as  j  -*■  » 
u 

fir 

Corr(VJZ.  ,  ,  VJZ  )  =  cos  (w)v  •  (dto)  -*•  -1. 

l~  ‘  1  L  TT  J 

But  this  means  that  for  every  tg 

Ed|^  =  P(dP^  =  l)-*-l, 

or 

E(d|J  +  1  ^  -*•  0,  j  -*■  , 

0  0 

and  apply  the  lemma.  u 

To  prove  that  a  process  is  not  strictly  oscillatory  it  is  sufficient 
to  show  that  Ed^  does  not  increase  with  j  for  some  t. 
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Example  1 .  Let  e  be  a  zero  mean  continuous  random  variable.  Define 


Ut)  by 
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«»%  m<t  »  » 


,«> 


Then 


and  so 


Zt  =  < 


,  t  <  0 

e]  »  t  =  1 ,2 

0  ,  otherwise. 


VZ 

V2Z 

V3Z 


£reo 


e0“el 


0 

0 

0 


1  =  Ed^  >  Ed^  =  P(eq  and  have  the  same  sign) 


and  the  process  is  not  strictly  oscillatory. 


The  last  two  examples  show  that  we  can  alter  a  sequence  of  random 
variables  and  even  a  strictly  oscillatory  sequence  and  render  it  non- 
oscillatory  in  our  sense.  This  leads  to  the  observation  that  an 
oscillatory  process  in  the  sense  of  Priestley  (1981)  need  not  be  strictly 
oscillatory. 


Theorem  4.  Suppose  a  process  { Z^}  is  an  oscillatory  process  with  the 
representation  (2).  Then  it  need  not  be  strictly  oscillatory. 

Proof.  It  is  sufficient  to  construct  an  example.  Let  {e^}  be  a  second 
order  Gaussian  stationary  process  with  mean  zero  and  correlation  func¬ 
tion  P^.  Let  {ct)  be  a  sequence  of  nonnegative  constants.  Define 


t  =  0,±1 ,. . . 


Then  is  oscillatory  with  respect  to  {ctelu)t}.  Assume  p2  =  0  and 
consider  the  time  points  7,  8,  say.  Then 


It  should  be  noted  that  by  a  proper  choice  of  c^,  may  still  be 
strictly  oscillatory.  For  example,  we  can  choose  ct  to  be  a  slowly  time 
varying  sequence. 


Consider  the  rate  processes  in  the  above  figures.  It  is  seen  em¬ 


pirically  that  from  a  certain  point,  say  Ng,  we  have 

D^N)  <  D2(N)  <  D3(N)  <  ...  <  1 ,  V  N  >  Ng, 

where  Ng  is  quite  small,  e.g.  Ng  =  100.  This  is  rather  expected  since 
by  Lemma  1,  for  all  N 

Dj(N)  <  Dj+] (N)  +  ^  surely.  (7) 

Now  fix  Ng  and  consider  Dj(N),  N  >  Ng.  Then  by  (7),  as  Ng  -+  °°,  Dj(N) 
for  N  >  Ng  will  tend  to  increase  with  j.  More  precisely,  (7)  entails 
that 

c.  =  lim  sup  D-(N) 

J  N-*oo  J 

converges  monotonical ly  with  probability  one  as  j  to  a  limit,  c, 

say,  such  that  0  <  c  <  1.  This  convergence,  though  useful  in  spectrum 
analysis  as  is  shown  in  the  next  section,  ignores  the  oscillatory  be¬ 
havior  of  the  process  for  low  and  moderate  N.  We  therefore  take  a 
different  route. 

Inspired  by  the  figures,  we  shall  assume  that  a  limiting  rdte 
process  exists  and  let  D(-)  be  the  corresponding  point  process. 

A  convenient  way  to  describe  the  monotone  convergence  of  the  rate 
processes,  is  to  prove  the  convergence  of  the  corresponding  point 
processes  D-|(*)»  D2(’ ),  ...,  obtained  from  strictly  oscillatory  pro¬ 
cesses,  to  D(*).  To  do  so,  we  follow  the  work  of  Kallenberg  (1976). 


'»*  '■*.  >i  >t  ki.  .i  *it,**|  "*«  ‘»r 


‘*1-  it. 


'il  'il.  il  'll  si 


ivwuiriivuvwiivTPvnivnivi  nr.nn»* 


For  p  e  N  we  write  pf  for  the  integral 


f(s)p{ds). 


For  each  fixed  f,  this  integral  is  a  mapping  from  W  into  R.  Let  f  be 
an  arbitrary  continuous  function  f:  S  -►  [0,°°)  with  compact  support. 
The  class  of  all  finite  intersections  of  N-sets  of  the  form 


{p:  s  <  jjf  <  t} 

where  s  and  t  are  real  numbers,  form  a  base  for  a  topology  on  N  to  which 
we  refer  as  the  vague  topology.  Then  Pn  tends  to  p  in  this  topology  iff 
bnf  -*•  pf  for  all  continuous  functions  f  from  S  to  [0 ,«>)  with  compact 
support.  Then  N  becomes  a  metric  space  and  we  can  think  of  the  D-(-) 
as  random  elements  in  (N, F).  By  definition  D-(*)  D(-)  if 

Eg(Dj)  -*■  Eg(D)  for  all  bounded  continuous  functions  g:  N  -*■  R. 

Observe  that  for  any  B  e  8  we  can  have  at  most  a  finite  number 
of  integer  points  and  thus  also  a  finite  number  of  zero-crossings. 
Therefore 

lim  lim  sup  P(D  ■  (B)  >  r)  =  0. 

p-K»  j~KX>  ^ 

Next,  by  definition  of  strictly  oscillatory  processes,  for  any  I  =  (a,b], 
a,b  finite,  we  have  by  monotonicity 


ED j ( 1 )  -  CD( I ) . 

Third  consider,  without  loss  of  generality,  the  higher  order  crossings 
°j,N  and  Dj+1  N*  Then  by  Lemma  1  »  there  exists  tQ,  Dj  and  Dj  +  1>  such 


-  »  h  •  _*  ^  .  I  .  \  .*■  .1  .  I  .  ,  U  i  . 


such 


and 
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D:  u  =  * 

JiN  j  tQ 

D.  =  n'  +  d(j+l) 
J+l.N  Dj+1  dtQ 


Dj  <  Dj+1  surely. 


By  strict  oscillation 


P(dp)  >  1)  <  P(dp'+1  ^  >  1), 


Therefore,  by  Lemma  1, 


P(Dj  +  dp)  >  1)  <  P(Dj+1  +  dp+1 )  >  1) 


or 


P(Dj>N.0)  >P(Dj+liN=0) 

and  so  the  sequence  P(D.  N  =  0)  converges  monotonical ly  as  j  for 
any  N.  By  appealing  to  Theorem  4.7  in  Kallenberg  (1976)  it  follows  that 
D  •  ( • )  converges  in  distribution  to  the  point  process  D(-)  corresponding 

J 

to  the  limiting  rate  process. 


Theorem  5.  Let  { Z t }  be  strictly  oscillatory  and  let  D(-)  be  the  point 
process  corresponding  to  the  limiting  rate  procsss.  Then  D-(*)  +  D(* 

J 

with  respect  to  the  vague  topology  on  N. 
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l 

' 

6.  The  Stationary  Gaussian  Case. 

Theorem  6.  Let  {Zt)“_i  be  a  zero  mean  stationary  Gaussian  process,  and 
suppose  to*  <  i t  is  the  highest  frequency  in  the  spectral  support.  If 
lim  D - ( N )  <  io*/tt  with  probability  one,  j  =  1,2,...,  then 

_  ★ 

limlimD.(N)=^—  a.s. 

J-+O0  N-k»  ^ 


Proof.  Let  v.  be  as  in  Theorem  3.  Then  from  Kedem  and  Slud  (1982) 

J 

v.  =>  -J-6  *  +  i<$  *  ,  j  °°.  But  the  Gaussian  assumption  implies  that 

j  2  -to  2  to 


cos  (ttED  .  (N) ) 
J 


rTT 

cos  (to)v  •  (dto) 

-TT  J 


COS  (to* )  , 


j  °°. 


Therefore  ED.(N)  to  /it  ,  j  -*•  °°,  uniformly  in  N. 

J 

By  strict  stationarity  we  can  define  the  a.s.  limit 


c,  =  lim  D,(N). 

J  J 

Then  by  Lemma  1  c.  <  c.+,  a.s.,  and  so  lim  c  -  =  c  <  to*/ir  a.s.  Thus 

-1  J  j-MJO  ^ 

by  bounded  convergence  E(c-)  -*  E(c),  j  -*■  °°.  Strict  stationarity  and 

J 

Fatou's  Lemma  yield 

E(c  • )  >  lim  sup  ED  .  (N) . 

J  j 

Since  Ec,  -*■  Ec  and  ED.(N)  to*/t r  , 

J  J 

E(c)  >  to*/iT  . 


But  c  <  (o*/tt  .  Therefore  c  =  <o*/tt  a.s. 


Coronary  2.  If  u>  =  tt,  then 

1 im  Tim  D  .  (N)  =  1 

j-K»  ff+OO  J 

Proof.  In  this  case  0  <  D.(N)  <  1  surely.  □ 

J 

Let  p^(l)  be  the  first  correlation  in  {V^}. 

00 

Theorem  7.  Let  { Z t } -j  be  an  ergodic  zero  mean  stationary  Gaussian  process 
Then 

limD.(N)  =  —  cos~l  (p^+^ (1 ) )  a.s. 

N-x»  J  71 

Proof.  The  proof  follows  from  the  fact  that 

P^  (1 )  =  cos(ttED1  (N) ).  a 

Theorem  7  shows  that  in  some  cases  we  can  expect  the  rate  processes 
to  converge  into  straight  lines,  a  fact  that  can  be  used  in  classifica¬ 
tion  and  discrimination  of  stationary  processes. 
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